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Abstract. The closest point of a linear variety to an external point is found by 
using the equality case of an Ostrowski's type inequality. This point is given in closed 
form as the quotient of a (formal) and a (scalar) Gram determinant. Then, the best 
approximation pair of points onto two linear varieties is given, as well as characterization 
of this pair of best approximation points. 

1 Introduction 

In this paper, we answer an implicit open question by Ky Fan and John Todd [31 
page 63]. We give a determinantal formula for the point where the inequality of 
the above referred to authors turns into equality, thusly obtaining the point of least 
norm of the intersection of certain hyperplanes. We present a result, in terms of 
Gram determinants, for the minimum distance from a certain linear variety to the 
origin of coordinates [Proposition 12. 1| . We note that this formula generalizes the 
one Mitrinovic [2llH] has given in the case of two equations. This best approximation 
problem was dealt with in [TD] , where the centre of (degenerate) hyperquadrics plays 
a decisive role. In [lOj . no answer in closed form was given. 

In this paper, we give a new proof of Beesack's inequality (p^. Theorem 1]; [SJ 
Theorem 1.7]), by following arguments used in [31 page 63, Lemma]. 

The Beesack's formula [Theorem 13. Ij gives the point of a general linear variety 
closest to the origin of the coordinates. We extend the formula of Beesack [1] 
Theorem 1] in order to get the nearest point of a linear variety to an external point, 
in H". When extending Theorem l3.11 we obtain the projection of an external point 
onto a general linear variety [Proposition 14. 1| . This Proposition 14.11 is used for 
getting the best approximation points of two linear varieties [Proposition 15. 1| . Also 
a characterizion of the best approximation pair of two linear varieties is presented 
[Proposition 5.2]. 

Our context is the Euclidean space IR", endowed with the standard unit basis 

{^1,1^2, ■ • ■ , ~^n) 
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and the ordinary inner product 



it 



UiVi + U2V2 + • • • + UnVn 



where 



it = ailf 1 + 021^2 H + an~^n = (ai, 02, . . . , a„) , 



The Euchdean norm 



l^ll 



• ~? is used and the Gram determinant is 



0(^1, "^2, • ■ --.^r) = det 



1^2 • 1^2 



1^2 • 1^1 



~^2 • 



1 < r < n. 



1^2 ■ ■■ itr'ltr 

(1) 

l,f2,-.-,tr) - if 

are hnearly dependent. See, for example, 



,^^) > and 



It is weU known that G ^2, ■ • 
and only if the vectors ~]ti,J^2, • • ■ , 1^ 
% page 132]. 

Some abuse of notation, authorized by adequate isomorphisms, is to be declared, 
notably the identification of point, vector, ordered set, column-matrix. 

This paper is organized in seven sections. In Section 2, we present and prove 
a result. Proposition 12.11 which answers an open question of Fan and Todd and 
make a remark concerning a formula of Mitrinovic. Section 3 is dedicated to a 
generalization of Proposition 12.11 this meaning that we study the projection of the 
origin onto a general linear variety. In Section 4, we deal with the projection of an 
external point onto a general linear variety. In the Section 5, we treat the distance 
between two disjoint linear varieties. We get and characterize the best two points, 
one on each linear variety, that are the extremities of the straight line segment that 
materializes the distance between the two linear varieties. An illustrative numerical 
example is presented in Section 6. Finally, in Section 7, we draw some conclusions. 



2 The minimum norm vector of a certain linear 
variety 

In this section, we state the Proposition 2.1, which solves an old open question of 
Fan and Todd. The proof makes use of a result of the mentioned authors. 

The next result [3, page 63, Lemma] gives the radius of the sphere tangent to a 
certain linear variety, as the quotient of two Gram determinants. 

Theorem 2.1 Lei "^1 , 1^2, . . . , it m be m linearly independent vectors in IR", 2 < 
m < n. If a vector ~af G H" varies under the conditions 

it i • Iz?' = 0, with \ < i < m — 1 



then _^ _^ 

. G'(ai, a2,..., a 

Lr {a I, a 2, ■ ■ ■ , a m~l, O, m) 

Furthermore, the minimum value is obtained if and only if lit is a linear combination 
of a I, a 2, ... , Om- 

For the sake of completeness and for later use in the proof of our Proposition 
2.1., we present here, essentially, the proof given by Fan and Todd [3l page 63, 
Lemma] . 
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Proof: For the vector 1^ satisfying conditions ©J we have 
G "a^2, • ■ • , "c^m, a?) = -G {1^1,1^2, ■ • ■ , 'f^m-i) + ( •1>!)G ("f?"i, 1^2 
Hence 



^m) > 0. 



-A- — 5> G f a 1, 



By hypothesis, the vectors 



1^ •l^ — 



Q, are hnearly independent, so 

G ("(^1, 7^2, ■ • • ,~<im,^) = 

It follows that 

G {it I, 1^2, ■ ■ ■ , ~a^m-i) 



if and only if "af is a linear combination of ~(t i, 1^21 ■ ■ ■ ,~(t 



GC^i,"?2,... 

, d fn—li d rn) 



(4) 



if and only if the vector is of the form 1^ = ailt i + a2~(t 



Now we are in a position for stating the equality case. A determinantal formula 
for the closest vector to the origin lying in a certain linear variety is given. 

Proposition 2.1 1. Let iti, it 2, ■ • ■ , it m be linearly independent vectors in M". 
The minimum Euclidean norm vector in M" satisfying the equations 

"^1 = 

"^2»^ = 

: (5) 





it m 
^ m 




= 
1 




1 by 










it I • it I 


it I • it 2 




itl • itm-l 




itjn-l • itl 


ltm-1 • it 2 
lt2 




a m — 1 * ^ 7n—l 

a m-1 


a m— 1 * a 

a m 


it I • it I 
it 2 • itl 


it I • it 2 
it 2 • "^"2 




a 1 • a ,„_i 
a 2 • a m-i 


• itjn 

■q^2 • "a^m 


itjn-l • itl 
itm • itl 


it m-l • it 2 
itm • it 2 




a m — 1 * a — 1 
a ?n * a 


^ m— 1 * ^ rn 



(6) 

where the determinant in the numerator is to be expanded by the last row, in 
order to yield a linear combination of the vectors iti, it 2, . ■ . , itm- 



2. Furthermore, 



—j, _ G ("o^i, "^2, • ■ • , "o^m-l) 
S • S — „ ._L =5 =5 ^ 



G(ai, a2,-.-, a m-i, a,„) 



Proof: Part 1. We look for the scalars ai,a2, ■ ■ ■ ,Q.m, such that the vector 1^ 
f amltm satisfies the conditions ([5]). 



ai a'l + a2'c?"2 + 
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For that end, we solve the system 
it I • it I ~^ 

1^2 • 'ctl 



it I • it 2 
~^2 • 



~t \ • ~t r. 
~^2 • ~^r. 



it I • itr, 
~^2 • "^^r; 









"0" 




a2 









OLm-\ 













1 



As the vectors "^2, • • • , "^m are, by hypothesis, hnearly independent, the 
determinant of the matrix of the above system, which is the Gram determinant 

G {it I, it 2, ■ ■ ■ , ~tm-l, ~<t m) , 

is non null. 

So, by the Cramer's Rule, we have 



a,; 



it I • it I 

~^2»~^1 


it I • ~ti-i 
'^2 • 







it I • 

"^2 • 


"^2 • 


m 
m 


ltm-l*~<tl ■■■ 
ltm»~tl 


'tm-l • '(ti-1 
itm • ~ti-l 




1 


~<t m-1 • "^i+l 
"^m • ~^ i+1 


m — 1 * ~t m 
~t m * 0^ m 




it I • it I 
"^2 • "^1 


it I • it 2 
~t2»~^2 




~^2 • "^m-l 










• ~^2 




~tm-l • ~tm-l 
* ~^m— 1 







with i = 1, . . . m. 

Here, for brevity, we introduce some notations: 

_ _ Gj^-^ 

(^i — Q 1 Oil a i — Q ' — Q ' 
G = G(lLi,~^2-,---,~tm-l-,~t ^ 



where 



Gi 



it I • it I 
it 2 • ~ti 

~t m-l • 



'S'l • iti^i If 1 • o;j+i 
7^2 • ~<ti-i 'ai'2 • "t^i+i 

m— 1 * ^ 2+1 



m } 1 



and the symbolic determinant 



7^1 • ~t \ 

~t2*~^\ 

~t m-\ • ~t\ 



~t\*~ti^\ • "^i+i 

~t2*~ti-\ 'S'2»'^i+1 

~t m-\ • "o^i-l ltm-1 • 



"5^2 • ~(tm 



"^1 • 'tm 
~t2 • ~tm 



We get, using these notations and rearranging in a suitable manner the terms 
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of the determinants, 



1^1 • it I 
it 2 • ~<ti 



it 2 • '^2 



a 2 ' a 



^ m — 1 * ^771—1 
a m-1 



7^ 



1 • a ,i 



a m — 1 * a 771 

a 7n 



~^2 • "c?"! 



a 1 • a 2 

"0^2 • 'S"2 



a 1 • a „ 

~^2 • "o^ri 



"(J'ni-i • ll 

a m • a 1 



"(^ „i • 1^2 



It 



1 • ~f^7n 

a • a 



Part 2. It is just sufficient to use (H)), in order to obtain || s' 



For computational purposes, we notice that, in the numerator of ©, the coeffi- 
cients of the vectors it i, it 2, ■ ■ ■ , "^777 are the co-factors of the elements in the last 
row of the matrix 



it I • it I 
it 2 • "^^1 



"S" 7(7-1 • ~<t I 



it I • it 2 
it 2 • it 2 



~<t m-1 • ~C?"2 



it I • it m-1 
it 2 • it m-1 



~<t m-1 • ~<t m-1 



it I • "5^,; 
it 2 • "S'r; 



a 777— 1 



• itr 



Remark 2.1 The particular case of Mitrinovic 

The determinantal formula given in ^1 of Proposition 2.1, for the least norm 
vector of the given linear variety is a generalization of the formula of Mitrinovic |^ 
page 67] and f% page 93] 

p p 

h J2 «f - Ofc J2 ^i^i 

(s-O im-Kk"") 

where (ai, 02, ... , Op) and (61, 62, ... , &p) are two non proportional sequences of real 
numbers satisfying 

p 

aiXi = and 

7=1 

3 The minimum norm vector of a general linear 
variety 

Here we treat the projection of the origin of the coordinates onto a general linear 
variety, so extending Proposition 12.11 The point where the sphere centered the 
origin is tangent to any linear variety is given in closed form by next relation ([5]). 
This result has been obtained in a different form and by another approach in py. 
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Theorem 3.1 ([I]) Let 

~^ 2t ■ ■ ■ m be linearly independent vectors in IR", 
with m > 2. The minimum Euclidean norm vector in M" satisfying the equations 

'^i • "a^ = Ci 

: (7) 

^ rn — 1 * — ^m — l 

with, at least, one non zero Ci^ i=l,. . . ,111, is given by the relation 



7^ 



whe 



7,. 



it. 



1 , . . . , m — 1 , 



a' r. 



Furthermore, 



G f a'l, a'2, 



, 0, tn — 1 



G [a' I, a'2, 



, a' 



-1, a 











a' 1 • a' I 


a'l • a'2 


a' 1 • a'm-i 


a'l • 0! rn 










a'm-i • a'l 


a'm-i •a'2 


a'„j_i • a',„_i 


^ m — 1 * <^ m 

— > 


7. 


^2 . 


a ra-\ 












a' 1 • a' 1 
c?2 • a'l 


a^i • 0^2 


a'i»a',„_i 

a 2 • 0! ra-\ 


o^i • a^m 


a' 2 • a' 2 


a'2 • a'™ 










a'm-i • a' 1 
^ — > 


a^i_i»a'2 • 


a'jm-\ •_^' m-1 


^ m— 1 * m 




a',„»a'2 


a • a 1 


a 9 a 



, (8) 



Proof: Performing elementary matrix operations, we turn into the form [ 
the last column of the augmented matrix of the system ([7]) 



Oil 

a2i 



ai2 

022 



ai„ 
o-2„ 



ai„ 
a2„ 



On 



Cl 
C2 

^m— 1 
Cm 



where 1?",; = (a^^jOi^, . . . ,ai„_i,ai„). 



(9) 
(10) 

(11) 

•• 1]^ 



4 Projection of a point onto a linear variety 

For dealing with this problem by taking into account the result of the preceding 
section, we use the fact that Euclidean distance is preserved under translations. 
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We are given a linear variety V and an external point Q. We perform a trans- 
lation towards the origin O of the coordinates: the pair {Q, V) turns into the pair 
{0,V'). We, then, apply Theorem 13.11 to the pair {0,V'). Finally, we undo the 
performed translation: we go back from the origin O to the point Q. We state the 
following 

Proposition 4.1 Let iti, it 2, ■ • • , l^m be linearly independent vectors in ]R", with 
m > 2. Then: 

1. The projection S of the external point Q := if = (qi, (72, • • • , 9ti) onto the 
linear variety V defined by 



it I • = ci 

1^2 • a? = C2 



(12) 



with, at least, one non zero Ci, i — \, . . . ,m, is given by 



s 



(13) 



where 



7' 



with 



^1.^2 



a m-1 • a 1 a „i-i •12 



7', 



"^=? — 

a"i . a" I 



a"i . a" 2 

^2-^2 



7^' mTt' 7^' mTt' 

a m-1 • a 1 a m~i ' a 2 
a m* a, I a m * 0, 2 



7' I •a/' m-1 



7^' mTt' 7^' mT^' 

a m-1 • a m-1 a m-l • a r: 

<^ ?n — 1 w m 



a"r. 



-lir 



a" I ma' m-1 
a'2 • a 'm-1 



a" I . a'\ 
7' 2*7', 



7^' mTt' m~^' 

a m-1 • a m-1 a m-i • a r 

m ^ ^ m—1 ^ m ^ ^ m 



i — I, . . . ,m — 1, 



(14) 
(15) 



and 



(16) 



2. For the distance, we have 
d\Q,V)^d\0,V') 



7' 



r(7 7 7 \ 

u \a I, a 2,---,a m-i I 
G (7 1,7 2, ■■■,7 m-1, 7 „ij 



(17) 



Proof: 
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1. We perforin a translatfon towards the origin of the coordinates, of the pair 
(Q, V) in order to get the pair (O, V). We have 



Replacing, in equations p^ . 1^ with x' +~<t, we get 



a 1 • a; 

~(^2 • a;' 



(1 jyi 9 X , 



(18) 



with, at least, one non-zero and = q — at • 

Now, by using relations ([8]), (|9]), ([TO)) and (fTTj) . we obtain the relations ([14 
(USD and (USD. 

Finally, undoing the translation, we have 
2. The Euclidean distance is translation invariant: 



d^{Q,V) = d\Q,S) = d\0,V') = 



-^2 (_r I a 1, a 2, • • ■ , a m-i I 



G l^a 'i, a '2, • ■ • , a 'm-i,a 'm^ 



5 Distance between two linear varieties 

In this section we deal with the interesting problem of finding the best approxi- 
mation pair of points of two given disjoint and non-parallel linear varieties Vi and 
V2. In other words, we are looking for the point on the linear variety Vi and 
the point S'2 on the linear variety V2 such that the vector S1S2 is, to within a sig- 
nal, the shortest one linking the referred to linear varieties. Here the main tool is 
the Proposition 133] This result is applied twice, just bearing in mind that, in the 
present case, the external point is either the generic point Gvi '■= gvi of the linear 
variety Vi or the generic point Gy^ := of the linear variety V2. 

Some notation is in order, for the sake of simplicity of the statement of our next 
result. 

We write the vector / e IR" the following manner: 

7 = (A, /2, . . . , A, h+iJk+2, . . . , /„) (/i, /2, . . . , A, 1^) e IR'' X iR"-'\ 

We state the main result of this paper 

Proposition 5.1 Let us consider two disjoint and non-parallel linear varieties Vi 
and V2 given, respectively, by 
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"^1 • "af = Ci 

"^2 • = C2 



14 := <^ 



(19) 



^mi — 1 
' mi ^ — "--mi J 

where l^i, '^2, • • • , "a^rm are linearly independent vectors in IR" and with, at least, 
one non zero scalar Ci, i = 1, . . . , mi, mi > 2, 



6 2 • = ^2 



^2 := 



(20) 



m2 — 1 * — ^7712 — 1 

^ m2 * ~ dm2 ; 



where ~^2, ■ ■ ■ , ~^ m2 o.'"'^ linearly independent vectors in IR" and with, at least, 
one non zero scalar di, i = 1, . . . , m2, m2 > 2. 

Let us denote 1^ = (xi, . . . , , Xm^+i, . . . , x„) G Vi as 1?" = ^xi, . . . , x^i , ^ ^ G 

) e V2 as V = (2/1, . . . e 



■ ; 2/m2 ) 2/m2 + l ) • ■ • ) 2/ra 



IR^ixlR"-"^ and if = (2/1,. 

Lei MS denote by [515*2] the shortest straight line segment connecting the two 
linear varieties Vi and V2 ■ 
Then 

1. The points Si G Vi and S2 € V2 are obtained through the unique solution of 
the overdetermined consistent system of linear algebraic equations 



(21) 



51 (if) = Gv, [ e' 

52 (t) = Gv, it) 



where: 



(i) Gvi := Gv^ {xm^+i,Xmi+2, ■■■,Xn) andGv2 i~ff) '■= (yms+i, ym2+2, ■■■,yn) 
are the generic points of, respectively, the linear varieties Vi and V2; 

(a) Si {If) := 5i {ym2+i,ym2+2,---,yn) 
and 

52 C ^ := 52 {Xmi + l, Xmi+2, Xn) 

are given, respectively, by 

Si{lf) = S'^{t) + Gv2{t) 
and 

S2(t)=S^(i)+Gv, (t); 
and where: 
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(iii) Si(lf) is given by 



7' I •J' I 



7 mi-l •a/' I 

^1 



a" I • a" I 



mx-\ •~a!' \ 



^ m 1 — 1 • <^ m 1 - 

^ mi — 1 



7'i • O^mi 



^ m 1 — 1 • rni 

^ mi 



fli'imi-i ai»a 



2*7' mi-1 



a 2 • fl' 



?ni 



a' 



-1 • mi-l O^r 



being 



' mi — 1 * ^ mi — 1 ^ mi — 1 * mi 
^ in\ • ^ mi— 1 mi • ^ ra\ 

(22) 



a' ai a z = 1, . . . , mi - 1, 

Wni 



^ rni ' 



mi 



mi 



m^/ij at least 
and 



one non zero c- — ci — • g\^, z = 1, . . . , uii, 



6^1.6^1 



7\ 



6"i.&"i 



7 



h" 

U r 



m2-l • b"l 



b"i 



6" 1.6" 



m2 



^"m2-l • &"m2-l 
m2 — 1 



&"i • b' 



1^' 



m2 — 1 

m2 — 1 



• &^m2 



1^1 



m.2 — 1 m2 — 1 



.6^ 



m2 — 1 



being 



m2 — 1 
6^m2 • 

(23) 



b I = b,~ -f- b r. 
^ 1 

^ m2 ~ "T/ b ra2 

"'m2 



1, . . . ,r7i2 - 1, 



wif/i, at least, one non zero d[ — di — b i • gvi , i = 1, . . . , m2 . 

2. The distance d(Vi, V2) between the two linear varieties is given by 

d{VuV2) 



S1S2 



Proof: Essentially the proof consists on dealing once at a time with the two linear 
varieties Vi and V2: 



1. finding the generic point of each linear variety; 

2. applying the Proposition HTT] 
In the following way: 
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(i) The generic points 

From the underdetermined system (IT^ . we can, without loss of generality, 
assume that the generic point Gy^ ■— 5Vi depends on the n— toi + 1 parameters 

We write 



Gvi = Gvi 



'^l('^mi+l; • ■ • ; '^n) 

('^mi + 1 1 • • ■ ; ) 
*^mi +1 



5^1- 



(24) 



Similarly, we write for the generic point 9V2 of the linear variety V2: 

yi{Vm2 + l, ■■■,yn) 



— Gv2 {~if) 



2/m2 (j/m2 + l 7 • • ■ 7 2/ri) 
2/m2 + l 



.9V2- 



(25) 



(ii) The application of the Proposition 14.11 

(a) Concerning the pair {Gv2, Vi), we get 



Siiym2 + l,ym2+2, ■ ■ ■ ,yn) — 



— ^ — ;> 
a mi-i • a 1 

^1 



a" I . a"i 



a mi-i • a 1 
a 'mi • a '1 



1 • O^mi-l 



^ 771 1 — 1 • m 1 — 1 

a 7711 — 1 



1 • 7771 



^ — > 

777 1 — 1 • 777 1 

7771 



a'l • a ',,11-1 
a '2 • a '7771-1 



a 1 • a 

a '2 • a '7711 



— 



— !• 



a 7771 -1 • a 7771-1 a „jj_i»a 



a 7711 'a 7711 — 1 



where 



a i 



It,. 



, mi - 1, 



a r 



with, at least, one non zero = q — a • gvj, « = 1, 



7 mi- 



ll 



(b) Concerning the pair {Gy^^V2), we get 



Pi •Pi 



Pm,-l»Pl 



P. 



Pi •Pi 
P2*Pl 



Pm,-l»Pl 
Pm, • Pi 



Pl^Pm,-l Pl'P'm, 

P mP p »p 

^ 1712 — 1 ^ 7712 — 1 ^ 7712—1 ^ 7772 
P 



7772-1 



6^1-&j,772-l 
&^2»&^7772-l 



6" 2 • b" 



7772 
7712 



^"m2-l • ^"7772-1 

P mP 

^ 7772 ^ 7772 — 1 



/,// - Ul 

^ 7772 — 1 7772 

P •P 

^ 7772 ^ 7772 



where 



1^1 . ^ 

i= hi- ^ = l,...,m2 - 1, 

1 

^ 7772 ~ ~Ti b m2 

with, at least, one non zero d^ = di — • g\^, i = 1,. .., m2- 

Essentially, the points S'l and 5*2 resulted from translations of the pairs 
(Gy^, Vi) and (Gvi,V2)- Undoing the translations, follows 

Si = S'l + 

5*2 = 5*2 + Gvi ■ 



We must get the unique solution of the ovcrdctcrmincd system 

'S'l(2/?772 + l) 2/7772+2) • • • ^ Un) = (x^i +1 , 2^7771 +2 , • • • iXn) 

'^2(^7771+1 ) *^777l+2) • • • ) ^77 ) = GV2 (j/7772 + 1) 2/7772+2, ■ ■ ■ , Vn) 

of 2n equations and the (n — mi + 1) + (n — m2 + 1) indeterminates 

^^7771+1? •^7771+2? • • • 5 ^77 5 2/7772 + 1? 2/7772+2) • • • ) 2/77* 

This system is consistent and has the unique solution 

1^7771 + 1 ' ^7771+21 • • • ) ^77' 2/7772 + 1 ' 2/7772 +2 ' ' " ' ' ^n/' 



(26) 



Hence we obtain 



Si=G 



•^1 (•^7771 + 1 ) ■ • ■ ' -^n) 

2;7ni (a^„ij+i, • • • , a;„) 













T* 












•^77 



and 



52 = G*y^ = Gv2 it*) 



2/l(2/m2+l'---'2'n) 

2/7772 (2/m2 + l' • • • ' 2/77,) 

2/7772+1 

2/;; 



2/7773+1 
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The assertion on consistence of system and uniqueness of the solution 
of system ((26)) is supported by results on existence and uniqueness of best 
approximation problems ^ page 64, Theorem 1] (4j page 45, Theoreme 2.2.5]. 



Remark 5.1 Some attention must be paid to the formulas \2°^) and \2S\) . In fact, 
we have 



s 



(27) 



i=i 



where A, Ai, i — 1, . 
and 



. , nil ttfe higher-degree polynomials in several variables j/m^+i, j/„i2+2, 



(28) 



where B, Bj, j = l,...,?7i2 are higher-degree polynomials in several variables 
However, from |^7| ) and i28\) we have 



7i = 7i{t)^J2^^^ (^)^ 



(29) 



where Lu 



, n, are first degree polynomials in the variables ym2+ii J/m2+2, ■ ■ ■ ,yn 



and 



(30) 



where L2i, i — 1, . . . , n, are first degree polynomials in the variables Xmi+i, Xm2+2, ■ ■ ■ ■, 
This question is worth a longer explanation. As follows: 

By performing the mentioned convenient translations on the systems 11 9\) and 
h2(]\) . we obtain two systems where the right hand sides are vectors whose entries 
are linear expressions in the parameters that are coordinates of the vectors Gvi and 
Gv2- By using arguments involving the uniqueness of (least squares) solution of 
a linear system by using the Moore-Penrose inverse, we assert that the solutions 
of the afore referred to systems are given in terms of such parameters. The best 
solution in the least squares sense of the system Ali — b is given |3 page 4-39] 
by = A^ b , where A^ stands for the Moore-Penrose inverse of matrix A. In our 
case, is a constant matrix, so depends on the parameters in vector h . 

Hence, 

Lliym2 + l,ym2+2, • • • , Vn) 
_^ L2{ym2 + liym2+2, ■ ■ ■ ,yn) 



Ln{ym2 + l,ym2+2, • ■ • , 2/n) 



-^2(^mi + l5 ^mi+25 ■ ■ ■ i ^n) 
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For the sake of clarity, we synthesize: 



Scholium Regarding the given linear varieties and without loss of generality, we 
can write 

■^1 (*^Tni +1 ; • ■ • ; '^n) 



*^mi +1 



Vo = < 



m2 + li 



,yn) 



2/m2(y m-2 + 1 7 • ■ • 5 2/n j 
2/m2 + l 



TJ12 + 1 : 



Hence we may write 



and 



Si — Gy^ — 



'-'2 — LriA — 



^mi + 1 



2/l(2/m2 + l.---'J/n) 
ym2 (ym2 + l' • ■ • ' 

ym2+i 



, . . . , y„) € H" 












*^7ni 




^r7ii+i 








2/1 




2/7712 




2/m2 + l 

. 2/ ;; . 



where 



V^^mi + li 2;mj4_2, ■ • ■ , X^, ym^ + i, 2/7712+2 

is i/ie unique solution of the overdetermined system 



>2/7l) 



A classical projection theorem [51 page 64, Theorem 1] g] page 45, Theoreme 
2.2.5] [21 page 64, Exercise 2] concerning the case of a point and a linear variety, 
leads us to a result on the projection vector connecting two linear varieties. It is a 
characterization of the pair of best approximation points, that may be useful when 
testing the accuracy of numerical examples. 

Proposition 5.2 Let Vi and V2 be two non-parallel linear varieties: Vi = Pi + Mi 
and V2 = P2 + M2, where Mi and M2 are subspaces of H" and Pi and P2 are 
fixed points in H". Then, the unique points Si € Vi and S2 G V2 form a best 
approximation pair {Si,S2) of the linear varieties Vi and V2 if and only if the 
two vectors whose extremities are Si and S2 are orthogonal simultaneously to the 
subspaces Mi and M2 ■ 
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Proof: We need just two facts: the definition of a vector orthogonal to a set of 
K" where a vector is said to be orthogonal to set if it is orthogonal to each vector 
of the set; and a projection theorem, where it is stated that the projection vector 
is orthogonal to the unique subspace associated to the given linear variety and not 
to the linear variety itself [5J page 64, Theorem 1] [U page 45, Theoreme 2.2.5] [U 
page 64, Exercise 2]. 
We have: 



1. S2 = ~S2 is the projection of Si := st onto the linear variety V2: hence S1S2 
is orthogonal to the subspace M2; 



2. 5i = it is the projection of 6*2 :— ~S2 onto the linear variety Vi: hence 5152 
is orthogonal to the subspace Mi. 



Notice that the vector S1S2 is not orthogonal either to the linear varieties Vi or 

V2. 

Finally, we have a result concerning the separating hyperplanes W) pages 105- 
106] and the smallest sphere tangent to the two linear varieties simultaneously. 

Corollary The smallest sphere S tangent to the linear varieties Vi and V2 is given 
by 



S= lit e H" 



^1 + ^2 

















1 




2 





and the supporting hyperplanes are 

H,^{lt e H" : (^1 - ^2) • - it) = 0} , j - 1, 2. 



6 Illustrative numerical example 

We are given two linear varieties. We exhibit the best two approximation points 
— one point on each linear variety — and show that the vector S1S2 is orthogonal 
to both the subspace Mi and the subspace M2 associated to the linear varieties Vi 
and V2, respectively, but not to the linear varieties themselves. 
Let the two linear varieties Vi and V2 be defined as follows 

Vi :— < overrightarrowai • "a?" = 1 (31) 
[ "^2»i^ = 2, 

with "^1 = (1, -1, -2, 1, 1) and ~^2 = (1, 1, -4, 1, 2); 

[ "^1.1^ = -10 
1/2 < ~t2»t ^ -20 (32) 
[ "^3-^-3, 

with 1^1 = (1, -1, -2, 1, 1), 1^2 - (-1, 1, -4, 1, 2) and 1^3 = (1, 1, -4,-1, 3). 
* 

Concerning the Proposition 5.1. 
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(I) The generic points Gvi and of the hnear varieties Vi and V2 are 

I + 3x3 — X4 



G 



Vi 



1^5 



2^3 
^5 



and 



23 
2 



-32/4 

hi' 
2/4 
2/5 



2 2/5 
- 52/5 
52/5 



(II) We perforin a translation along the vector Gv2 {yA^y^jO = O—Gv^ (2/4; 2/5); 
the linear variety ¥{ is obtained by replacing 1^ in the relation pip with 

X ■ 



G 



V2- 



Vi:= 




22/4 + 2/5 



getting 



-51(2/4,2/5) 



7 ^ 


^2/4- 


^2/5 


131 
21 


-§2/4- 


^-§2/5 


4 
21 


^§2/4- 


-§2/5 


15 , 
7 ^ 


^2/4- 


^2/5 


2 

21 


§2/4 + 


^2/5 



and 





- 191 

14 


+ 


§2/4 - 


§2/5 




221 

42 


+ 


63 2/4 


126 i/5 


Gv2 — 


101 

21 


+ 


§2/4 + 


126 y5 




15 
7 


+ 


§2/4- 


^2/5 




2 

21 




§2/4 + 


63 2/5 



Mutatis mutandis: 



(33) 



(III) We perform a translation along the vector Gvi{x3,Xi,x^)0 = O — 
Gvi{x2.,Xi,x<^); the linear variety is obtained by replacing in the 



— )• 



relation ([5^ with y + G 



^i-y[ = -ll 

:= <^ 1^2»?" = -19 + 6a:3-2x4 
b 3 • y' = 1 + 2x4 — 2^5 



3X5 



(34) 
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getting 



6'2(a;3,a;4,a;5) 



633 _ 366„ I 
209 209*^3 

35 I 12 _ 
19 ^ 19*^3 

674 _ 150 , 
209 209*^3 ^ 

_1371 , 2^ _ 
209 ^ 209-^3 

172 _ 42. , 
209 209*^3 T 



192 , 148 
209 4 209 5 

— g|a;5 

41 I 159 
■^4 T ^lo-tS 



209 



209*^4 418-^5 

70 7_ 

209 4 209 5 



and 



<S'2(a;3,a;4,a;5) = 82+ Gvi = 



418 



261 
209 



X3 



17 
209 



X4 



418*^5 



89 , 31 _ J_ _ 17 

3g -r iqJ-3 iqJ-4 iQ-i-s 



19 



19" 



19" 



674 , .59. , «_ , 159 
209 ^ 209*^3 ^ 209*^4 418-^5 



1371 
■ 209 



+ ii^3 + 



18 
209 



a;4 



418" 
418-^5 



172 _ 42_ , 70_ , 202 
209 209 3 ' dno"^* ^ 9nQ'*'o 



209" 



209" 



(IV) Solving the system 



'S'i(y4,y5) = Gvi(a;3,X4,X5) 



wp nhtain r* - — r* - -iI65 * _ 1489 * _ _3560 * _ 453 

we ODiam — — g^g , x^, — ^^a ' Vi — 509 ' ~ 212- 

Hence, using 



r 3 

2 



+ ix% 



x\ 



2-^5 



2-*'5 



X4 



and 



we, finally, obtain 



G 



V2 



¥ + ?4 



2 
23 
2 



5 + \yl 



2Vb 

I2/4 ~ 12/5 

1 * 
2^5 



yl 



77 

16 

57 
212 

837 

848 
4765 
848 
1489 
424 



and ^2 



55 

16 

469 
424 

3169 
848 
3560 
509 
453 
212 



The distance between the two varieties is given by 



d{VuV2) = \\SiS2\\ = 



2174 
559 ■ 
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Concerning the Proposition 5.2. 
Let us consider 



= Pi + Ml := 



1-3-1 

2 






■ 2>X3 - X4- |a;5 


> 


1 
2 






X3 - \xr, 







+ < 




X3 


: 0:3, 0^4, e IR > 









Xi 











Xr, 


> 



and 



V2 = P2 + M2 := 



- 23 - 

2 






23 
2 






5 


+ < 
















k - 



Vi - 

4 

32/4 

|y4 



2?/ 
1 



2?A=5 
5^5 



2/4 
2/5 



J/4,2/5 e IR 



ai) The vector S1S2 is orthogonal to the unique subspace Mi associated to 
the hnear variety Vi. Consider the arbitrarily fixed vector 

3, 



^1 



We have 



3X3 — X4 

X3 
X4 

X5 



e Ml. 



S1S2 • = 0. 



0:2) The vector S1S2 is orthogonal to the unique subspace M2 associated to 
the linear variety V2. Consider the arbitrarily fixed vector 

1, 



V2 



2/4 - 
4 

3^4 

hy4 



2-2/5 

- ^2/5 



2/4 
2/5 



G M2. 



We have 



^1^2 • ^ 



0. 



Pi) The vector S'iS'2 is not orthogonal to the linear variety Vi. 
Take the fixed vector 



e^i. 



We have ^1^2 • tli = -1.3750 ^ 0. 
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(32) The vector S1S2 is not orthogonal to the linear variety ¥2- 
Take the fixed vector 



^2 



23 
2 

71 

6 

24 
3 

1 

2 



e ^2. 



We have ^iS-s •U2 = 13.7500 7^ 



7 Conclusions 

In this paper we presented a determinantal formula for the point satisfying the 
equality condition in an inequality by Fan and Todd (we answered the implicit old 
open question in [21 page 63]: to get a closed form for the minimum norm vector 
of the given linear variety). In a previous paper [TU], we got, by using the center of 
convenient hyperquadrics, the point where the inequality ([3]) turns into the equality 
©. 

Here, we also restated a determinantal formula for the point of tangency between 
a sphere and any linear variety. 

Furthermore, we obtained the projection of an external point onto a linear va- 
riety as a quotient of two determinants. Subsequently and consequently this result 
was extended for getting the best approximation pair of two disjoint and non par- 
allel linear varieties. A characterization of this pair of best approximation points is 
offered. 
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